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1 Introduction 

We consider a new system of differential equations which is at the same time 
gradient and locally Hamiltonian. It is obtained by just replacing the factor 
i by i 2 = — 1 in the equations of interaction for N point vortices, and it 
is interpreted as an interaction of N point sources. Because of the local 
Hamiltonian structure and the symmetries it obeys, it does possess some of 
the first integrals that appear in the N vortex problem. We will show that 
binary collisions are easily blown up in this case since the equations of motion 
are of first order. We just pass to relative coordinates and then a convenient 
change of time scale. This method may be easily generalized to the blow up 
of higher order collisions. 

The author is indebted to Jesus Mucino for conjecturing that a sort of 
harmonic conjugate of the Hamiltonian function for point vortices could be 
a first integral of the vortex system. It turns out that this is not true, 
but nevertheless motivated the introduction of the model for a system of 
point sources. Furthermore, (M. Hampton, 2008) proposed to take complex 
vorticities instead of the standard real ones, which allows us to generalize the 
model to cover simultaneous interactions of sources and vortices. 
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2 Statement of the problem and main prop- 
erties 

The Hamiltonian for a system of planar point vortices is given by 

N 

h = - ^2 r fc r,in|^-^|. (i) 

k=i,i>k 

If we construct the gradient vector field corresponding to H, we get instead of 
the differential equations for N point vortices as in (Aref, 2007 and Newton, 
2001). 

N N -p 

E Zk — zi . sr^ J- 1 
Tl '\ 12 = 1 1^ ~ ~> 
\Zh — Zi\ z Z k — Zi 

Ijtk 1 1 l^k K 1 

those corresponding to the interaction of a system of planar point sources, 
i.e. 

N 7 7 N TV 

T k z k = - £ T^-^^- = - £ JiEL. (2) 

This velocity vector field X s can also be represented as a locally Hamiltonian 
vector field for the one-form 



N 



a 



= -d $]r fc r,axg(z fc -z,) , (3) 



J>fc 



i.e., «x s ^ = -a. 

This a is globally defined in C w \ A where A = {(zi,--- ,z^)\zj ^ 
ZkVj, k}, and in fact it is a closed but not exact one-form, since the function 
in between the parenthesis is only locally defined. This is because we can 
rewrite 

N 

a 

i>k 

and each polar angle 6 k i = aig{z k — zi) is defined only modulo 2n in C \ {0}, 
although its exterior differential d9 k i is globally defined there. 
One can easily verify that the linear momentum 



^r fc r,d(axg(z fc -z,)), (4) 



^r fc z fc , (5) 



k 

2 



is a first integral of system fl2]), so if J2 k Y k ^ we can define an equivalent 
center for the sources by zq = Z/^2 k Y k . On the other hand, the angular 
momentum A = Y2 k Y k z k x z k is an invariant of the system. In fact, we have 

Proposition. 1. The angular momentum A for the point sources system is 
identically zero. 

Proof. If • denotes the standard scalar product of R 2 identified with C, we 
can write 



^ Y k z k x i fc = Y k {iz k ) ■ z h = - r fc (iz fc ) ■ ( ^ 1^ ^ J = 




since Y k Yi(iz k ) ■ (z k - z t ) + Y k Yi(izi) ■ (z t - z k ) = 0, for any k^l. □ 

In some sense, the reason for this result is that we are dealing with a 
central velocity vector field. So, there is no simple way to introduce rotations 
into the system. 

Conjecture. . Along any fixed solution of the polar angles Q k \ undergo 
variations smaller than 2ir, so that g = — ^2^ >k Y k Yfi k i is a first integral of 
the system along the given solution. 

However, the moment of inertia / = J2k^k\zk\ 2 = ^2 k YkZkZk, which is 
a first integral in the case of iV-vortices, is not always a first integral for 
iV-point sources. Indeed, 

I = y Y k Z k h + ^ Tk ZkZk 
k k 

-p -p z k — Zi \ ^ _ \ ^ z k — Zi 

= -^r fc r z , 

since z k {z k — zi)+zi(zi — z^) = \z k — zi\ 2 . So, / is in general a nonzero constant. 
Hence, only the two point source problem is always integrable since it has 
enough first integrals. In order for the integrability to make sense, we have 
to definitely know that we are dealing with a globally Hamiltonian vector 
field. But this is the case for N = 2, since the angle 8u is not only bounded 
but constant, as it will become also clear from the analysis in the following 
section. 
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Proposition. 2. IfJ2i<k^k^i = 0, then I is a first integral. IfJ2i<k^k^i ^ 
0, then all solutions are gradient-like with respect to I. 

This implies that only when J2i<k ^k^i = is there any chance of having 
periodic solutions. 

Comparing with the iV-vortex problem, we recall that in that case if 
the Tfc denote vorticities, the angular momentum takes exactly the value 
J2i<k ^k^i, which is called the virial, while / is a first integral of the system. 

Also, collisions are rather scarce in the iV-vortex problem and depend on 
the values of the vorticities. In contrast, collisions are very likely to occur in 
positive or in negative time for a system of iV sources. There is no simple 
way to continue solutions beyond collisions, unless we agree that any time a 
collision occurs the sources involved merge into one source whose intensity 
is the sum of the colliding source intensities. As a result, the solutions will 
only be continuous but not smooth when collision occurs and the dimension 
of the phase space jumps correspondingly. 



3 Example: A single point source and two 
point sources. 

The differential equations for a single source are 

r r 

z = -—z = --, (6) 

\z\ 2 z 

where z£C \ {0}, or in terms of real and imaginary parts 

Tx . Ty 

x + y x + y 

where — T is the intensity of the point source. When — r < the vector field 
points radially towards the origin, so that it is actually a sink; when — r > 
it is really a source, since the fluid runs away from the origin. 

In this case we can take u = dx A dy and a = — T ~ y ^ 2 x ^ 2 dy . The one form 
a is closed but not exact on C \ {0}, since we can write a = — Td9, but 
only when the polar angle 9 is defined in an open interval of length smaller 
or equal than 2n. 
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The equations of 2 point sources of nonzero real intensities Ti and T 2 are 
given by 

• _ r Zl ~ Z2 ■ r z 2 — Zi /„\ 
\Z\ - Z 2 \ 2 \z 2 -z 1 \ 2 

In terms of relative coordinates z = Z\ — z 2 , they reduce to the single source 
equations of intensity — (Ti + T 2 ): 

i = -(r 1 + r 2 )^. (8) 

If — (Ti + T 2 ) < this is an equivalent sink, so that the 2 point sources 
approach each other. On the contrary, if — (I^ + T 2 ) > it is an equivalent 
source and the 2 point sources move away from each other. Finally, if Ti + 
T 2 = 0, then z is constant, so that the two sources will move with the same 
velocity and keeping a constant distance, as we can see from ([7]). 

In the cases where Ti+r^ ^ we verify that indeed / is not a first integral. 
All the solutions will then either end in collision or begin at collision in finite 
time. Indeed, equation ([H]) may be written as 

^ t \z\ 2 = zz + zz = -2(r\ + r 2 ), 

which implies \z\ 2 = — 2(Ti + T 2 )t. 



4 Regularization and blow up of binary col- 
lisions of point sources 

In this section we show how to desingularize any binary collision of point 
sources. 

We begin with the single point source which as we saw in the last sec- 
tion is equivalent to a two point source. Upon multiplication of the equation 
by z, we get zz = —T, suggesting the change of time scale — z, producing 
the new equation 

dz 
dr 

in terms of the time r. Moreover, we may assume zel, since any motion 
is radial. This equation can be integrated as z(t) = — T(t — r ) + zo, and 
from dt = z(r)dT we get the quadratic equation relating both times t = 
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-r^ + (Tr + zo) = -|r(r - r - z / r ) 2 + |r(r + ^ /r) 2 . This may be 
considered as a regularization since the collision is not a singularity any more 
and it is reached in finite time. However, the physical time is bounded and 
it is not defined beyond the collision time value t = ^r(r + z /T) 2 , which 
is reached when r = tq + zq/T. This is clear, since in contrast with the 
gravitational 2-body problem, we would not expect here the motion to be 
defined beyond a collapse of two sources since they merge into a single source. 
If we multiply instead equation (Q by \z\ 2 we get a new time rescaling 
= \z\ 2 which produces the differential equation 



ds 



— = -Tz 

ds 

whose integration gives an exponential solution z(s) = e~ Ts . This new desin- 
gularization corresponds to a blow up, since the singularity z = is now an 
equilibrium point and so it is reached in infinite time. 

Because for 3 or more sources collision does not occur in general along a 
straight line, it is the blow up which can be used for desingularizing binary 
collisions. We will illustrate the transformation for the case of 3 sources. 

The differential equations for 3 sources located at zi, z 2 , z 3 G C with 
respective intensities r 1; r 2 ,r 3 are given by 

r 1 



1 


Zl 




Z2 




1 




Z2 




Zl 




1 





3~ —■> 

Zl - Z 3 



Z2 = -Ti- — --r 3 - (9) 

z 2 - z 3 
1 

-r 



2~ Z - ■ 

Z3 — Zl Z3 — Z2 

Passing to the relative coordinates £ = z$ — zi, r\ — Z3 — Z2, we can write 
z 2 — Zi = £ — 77, getting the following system of differential equations for £, 77 

• = _ Ti + r 3 r 2 r 2 



. = _£i _ r 2 + r 3 _ i\ 
v C v fj-z' 
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|£| 2 yields the blow up system 



ds 



(r\ + r 3 )£ 





i\£- 



(r 2 + r 3 )|ei 2 r^i 2 



which is already regular at binary collision £ = 0. If we want the solution in 
terms of the original coordinates when r = I\ + r 2 + T 3 is nonzero we need 
to solve the system 



If we want to blow up at the same time binary collisions £ = and rj = 0, we 
change the rescaling by ^ = |£| 2 M 2 - This situation is generalized as follows. 

Proposition. 3. For the general case of N sources we change to N — 1 
relative coordinates £1,^2, • • • , £jv-i an d take a time scale factor j| which is 
the product of each of the |£j| 2 corresponding to the collisions £j = we want 
to desingularize. 

Passing to relative coordinates is equivalent to factoring out the transla- 
tion symmetry corresponding to the first integral Z. One can show that for 
N sources, only N — 1 suitable chosen relative coordinates generate all the 
remaining ones. 

If we try to formally apply this regularization and blow up transforma- 
tions to the differential equation z — for the velocity vector field of one 
vortex, we get in the first case zz = iT, but rescaling does not have any sense 
here if we want to keep a real time; but by taking real parts this equation 
yields jj\z\ 2 = 0. In the second case we get ^ = \z\ 2 z = iTz, which upon in- 
tegration gives z(s) = Ke tVs and finally j| = \z\ 2 = K 2 and so t — K 2 s + t . 
In both cases, we see that solutions are concentric circles as it is well known. 



T 1 z 1 + T 2 z 2 + T 3 z 3 = Tz 



-z 1 + z 3 = £ 

-Z 2 + Z 3 =Tj, 

where zq is the equivalent source center for the system, getting 



^i = ^o + (£ + Wr-£ 
z 2 = z + {t + v )/r - v 



zs = zo + (£ + ri)/T 
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5 Interaction of point sources and point vor- 
tices 



We can generalize the models for interaction of point vortices and of point 
sources to a model where we consider simultaneous interaction of point vor- 
tices and point sources. 

The idea is to formally write the differential equations for point vortices 
where we assume now that the vorticities are nonzero complex numbers, that 
is 



TV 

E F z k — z l 
Vl '\ 12: 
\Zk - zA 2 



l^k 

with Tfc G C \ 0. This new velocity vector field X can be written as the 
sum of a Hamiltonian vector field Xi with respect to the function H = 
— X)fcLi i>k R- e (r/c)Re(r;) In \z^ — zi\, plus a gradient vector field X 2 with re- 
spect to the function G = — J^fcLi i>k ^ m (Xk)^ m (Xi) In \ z k — %i\ if we decom- 
pose the system as 



N N 



Zk 



= * E - E Mro^' 
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The vector field X 2 is also locally Hamiltonian with respect to the one-form 
a = - J2ilk ^k^i d(arg(^ fc - z t )), globally defined in C N \ A. Here, Re(I\) 
represent the vorticities and — Im(rjt) represent the source intensities. 

If not all the real numbers Re(Ffc) are nonzero, we reduce correspondingly 
the dimension of phase space in order to define the corresponding symplectic 
structure to produce the Hamiltonian vector field X 1 . 

5.1 Example, 2 point source- vortices 

The equations of 2 point source- vortices of nonzero complex intensities r x 
and T 2 are given by 

Zi — z 2 . z 2 — z 1 

Zl = -r 2] TT, Z 2 = -IV 7Z. (10) 

\zi-z 2 \ 2 \z 2 -z 1 \ 2 

In terms of relative coordinates z = z 1 — z 2 , they reduce to the single source- 
vortex equations of intensity — (T 1 + T 2 ): 

i = -(r 1 + r 2 ) T 4 ? . (11) 



\z 
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The velocity vector field corresponding to these equations is considered in 
Chapter 11 of (Needham, 1997) in a different context. This time equa- 
tion ffTTl) implies \z\ 2 = — 2Re(Ti + T 2 )t + 7q. But in fact, if we write 
z = re tB then z — (r + ir6)e ld or zz = (rr + ir 2 6) = — (Fx + r 2 ), i.e. 
rr = — Re(ri + r 2 ),r 2 # = — Im(Ti + T 2 ). Integrating, we get 

r = ^/-2tRe(ri + r 2 )+rg 

and 

f -lm(T 1 + T 2 )dt 
J -2tRe(r 1 + T 2 ) + r 2 + "° 

This corresponds to a logarithmic spiral of the form r = Ke ad . So that 
a collision occurs in finite time with infinite spiralling, that is, while the 
argument 8 remains unbounded. If RefTx + T 2 ) = 0, then r = r is constant 
and 9 is linear in t, corresponding to an equivalent pure vortex. Likewise, 
if ImfT! + T 2 ) = 0, then 9 = 6 is constant and r 2 is linear in t, which is 
equivalent to a pure source. 

If we apply now a blow up to equation ffTTj) . we get \z\ 2 z = — (Ti + 
r 2 )z. The change of time scale ^ = \z\ 2 gives ^ = — (I^ + T 2 )z, which is 
integrated as z(s) = Ke'^ 1 ^ 3 = Ke- Rc( - Tl+T ^ s e- ilm( - ri+T ^ s , giving another 
description of the logarithmic spiral in terms of time s. We remark that 
an infinite spiraling while going to collision in finite time occurs also in the 
Manev problem which is an approximate model for the Schwarschild equation 
in general relativity. 
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